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The possibility of describing the behavior of a QQ pair in a hot plasma in terms of an effective 
potential is investigated. It is shown that as long as medium effects can be embodied in a gaussian 
action, like in the QED case, the QQ propagator obeys a closed temporal evolution equation 
whose large-time behavior is governed by an effective potential. The latter, beside screening, 
displays also an imaginary part related to collisions. 
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1. Introduction 

More than 20 years ago Matsui and Satz ^ proposed the J/y suppression due to Debye 
screening as a signature of the onset of deconfinement, triggering a lot of experimental efforts 
to measure its yields in high-energy nucleus-nucleus collisions. On the theory side in the last 
few years quite precise lattice data became available, both for the QQ free-energies [^] and the 
in-medium correlators of quarkonia [|0. This, hopefully, should make possible a more solid quan- 
titative study of the problem. Quite surprisingly the spectral functions extracted from the tem- 
poral correlators of quarkonia turned out to display well-defined peaks, reflecting the presence of 
bound/resonant states till temperatures of order 2T C , at least in the s-wave channels. Various at- 
tempts were done in order to interpret both the above large melting temperatures [Q and the 
charmonia correlators themselves [||, [7|] in terms of screened potential models. The latter repre- 
sents a quite economic way of accounting for medium effects, which is often employed to describe 
the interaction between charged particles in a polarizable medium. However, beside checking the 
numerical agreement of the findings obtained in the two approaches, it would be of interest to see 
such a medium-modified QQ potential arising from a first principle calculation. Different papers 
appeared quite recently addressing this issue [§,|,[T(j]. Here we wish to answer a few very general 
questions concerning the description of a QQ pair in a hot plasma: whether it is possible to give a 
solid theoretical basis to the concept of an effective in-medium potential; if so, which is its link with 
the QQ free-energy obtained from the imaginary-time propagator of a static pair of heavy quarks; 
finally, whether, beside the screening of the interaction it is possible, within the same framework, 
to account also for other effects, like the collisional damping. 



2. The QQ propagator in the complex time plane 

We start our investigation from the following propagator 

G > (^r 1 ;^r 2 |0,r;;0,r^) = ^(^r 2 )vA(^r 1 )vA t (0/Jz t (0^)) = (/(^r 1 ,r 2 )7 t (0,r / 1 ,r^)) (2.1) 

of a QQ pair created at time and annihilated at time t. From its spectral decomposition it follows 
that G > (t) is an analytic function of the (complex) time t in the strip — /3 < Im t < 0. In the case of 
static quarks the above propagator reduces to 

G^ =oo (f J ri;f l T2|0,ri;0y 2 ) = 5(n-r' 1 )5(r2-i / 2 )G(t,r 1 -r 2 ). (2.2) 

The analyticity of G > (t) allows to follow its evolution for imaginary times. In particular, in the 
static case, its value atT = j6=l/T 

G(t = -iP,n -r 2 ) = e -^QQ^-^) i (2.3) 

gives the change of free-energy occurring once a QQ pair is added into a finite-temperature medium. 
The latter is also the quantity provided by the lattice-QCD simulations of Polyakov line correlators 1 
[11], hence its interest in the present contest. 



'The evolution of a static quark from t = to T= j3 is in fact described by a Polyakov line. 



2 



In-medium quarkonium correlators 



Andrea Beraudo 



In [||] the authors proposed the following strategy to properly define a real-time potential. One 
assumes that, for static quarks, G > (t) obeys an equation like (with simplified notation) 

[id t -V(t,r)}G> =oo (t)=0. (2.4) 

One can then evaluate G > (t) at the lowest order of the HTL resummed perturbative expansion 
as G > (t) = G^(t) + G%Jt) + . . . , plug it into Eq. (|2T4|), identifying in such a way the leading 
contribution to the effective real-time QQ potential 

V(t,r)=V^(t,r) + ... (2.5) 

The latter is then employed in the dynamical problem with finite-mass quarks, 

[id t -T-V^(t,r)}G > (t)=0, (2.6) 

implicitly assuming that also in this case the pair propagator obeys a closed Schrodinger equation. 

The above procedure is far from trivial. In general in fact the temporal evolution of a n — 
particle propagator is part of an infinite hierarchy of coupled equations. Is it possible to find some 
physical systems which share many relevant features with the QGP, but for which one can get a 
closed Schrodinger equation for the heavy-pair propagator G > (t) from a first-principle calculation? 
A hot QED plasma of photons, electrons and positrons will answer to our purposes and will be the 
subject of the next section. 

3. The QQ propagator in a hot QED plasma 

We concentrate on the case of a static pair and in order to evaluate its evolution in the complex- 
time plane we proceed as follows. The propagator in a given background configuration of the gauge 
field is given by the product of two Wilson lines 

G A (f,ri,r 2 ) = exp (igj^ dt'AQ(r x ,t')\ exp (~igJ Q dt'A (r 2 ,t')\ = exp (ij d 4 zJ ll (z)A^{z) 

(3.1) 

One should then average over all the possible field configurations with an action accounting for 
medium effects: 

G(t,r l -r 2 )=Z- 1 J [@A]G A (t,r h r 2 )e iS[A] . (3.2) 

For the latter a convenient choice is the HTL effective action, which allows to properly include the 
most relevant medium corrections to the propagation of long wave-length modes (A so f t ~ 1 /gT). It 
can be expressed in terms of the time-ordered (along the usual Schwinger-Keldysh contour C in the 
complex-time plane described in ^) photon propagator 

D MV (x -y) = i d c (x° - y°) (A„ (x)A v (y)) + i d c (y° - x°) (A v (y)A^ (x)), (3.3) 

taken in the HTL approximation, and reads 

S H c TL [A] = -\ J c d 4 x j c d*yA»{x) (D- i fJ v L (x-y)A v (y). (3.4) 

Being the above action gaussian allows to perform the functional integral exactly getting 



G(t,ri -r 2 ) = exp 



l -f c d*x j c d*yJ»{x)Df v \x-y)J\y) 



(3-5) 
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3.1 Real-time propagation 

A QQ pair propagating along the real-time axis is described by the current 

J»( z ) = 8» e(z )e(t-z )[g8(z- ri )-g8(z-r 2 )], 

which leads to 



G(t, ri 



-exp 



d(0 
2k 



dq l-cos(fiw) f iq (r 



(2k 



(O- 



Doo((0,q) 



(3-6) 



(3.7) 



where the FT of the time-ordered propagator (which we take in Coulomb gauge) can be expressed 
in terms of the HTL photon spectral function as follows [^] : 



D 00 (w,q) = — + 



dq° p L (q°,q) 
2n q° - (co + iri) 



+ ip L ((0,q)N((o) 



(3.8) 



The above exact exponentiation, occurring in the case of a gaussian action, allows an immediate 
identification of the real-time potential V(t,r) defined implicitly in Eq. (2.4). It is of interest, in 
particular, to consider the large-time behavior of the above propagator, which results governed by 
the static limit of the HTL photon propagator. One has: 



lim [id t -V eS (r l -r 2 )}G(t,r l -r 2 )=0, 
with the in-medium effective potential we were looking for given by 



(3.9) 



V eff (r! 



dq 



1 _g''q-( r i- r 2) 



KmffT 



lq 2 + m 2 D 



|q|(q 2 + m2)2 



4k 



,-m.nr 



g 2 T 



(3.10) 



The above potential, arising here from a first-principle calculation, embodies medium effects both 
in its real part, with the screening of the interaction and the heavy-quark self-energy correction (its 
finite value requires subtracting the vacuum Coulomb self-interaction), and in its imaginary part 
related to the soft collisions with the plasma particles suffered by the heavy quarks. Note that the 
self-energy correction is crucial to ensure that medium effects vanish for very small separations 
and that an effective potential of the kind 



Veff(r) 



■ami) 



a 



-trior 



(3.11) 



is known in solid-state physics as Ecker-Weitzel potential [O], often employed in the study of 
excitons in semiconductors. 

In [Jsj] the propagator in Eq. (2T) was given a gauge invariant definition by joining the heavy 
quark fields with two Wilson lines. However for the study of the t — > oo behavior this is of no 
consequence since only the diagrams with all the vertices attached to the long sides of the resulting 
Wilson loop give the leading contribution. The situation is different in considering the imaginary- 
time evolution, which occurs only till T = j3. The latter will be addressed in the next section. 
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3.2 Imaginary-time propagation 



We now consider the QQ propagation along the imaginary-time axis. It can be expressed in 
terms of the Matsubara (HTL resummed) photon propagator 



G(-iT,ri-r 2 
obtaining 

G(-rr,ri-r 2 



exp 



g 2 dx'\ dx" 



dq 



l-e iq - {ri - r2 AA L {x' 



(3.12) 



x exp ig 



exp 

dq 



dq 



1- 



(2»)3 
_^'q-(ri-r 2 ; 



^•(n-r,)^ fzl + [ d i° Pdq°,<i 



2k 



dq p L (g°,q) 



2k (q°) 



,0\2 



(^ T -l)(l+A^°)) X 



x exp < g 



dq 



l_ £J ''q-(r 1 -r 2 ) 



dq PL (q°,q) 
271 (q ) 2 



\)N(q° 



Its value at T = j8 



G(-/j3,r 1 -r 2 ) = exp^- J 8g 



dq 

J2ny 



q 2 + m^ 



(3.13) 



(3.14) 



allows (after subtracting the divergent vacuum contribution to the heavy-quark self-energy) to iden- 
tify the QQ free-energy, getting: 



g e 



g 2 m D 



4k 4k r 



(3.15) 



Note that the correlation function from which the latter was obtained, and which represents a sim- 
plified version of the usual Polyakov-line correlator evaluated in lattice-QCD, is indeed a gauge- 
dependent quantity. This issue was carefully examined in [13, 14], where the usual strategy of ex- 
tracting an in-medium effective potential from the lattice color-singlet free-energy was criticized. 
What we can say after our investigation is that, at least in the case of a QED plasma for which an 
exact exponentiation holds, the QQ free-energy evaluated in the Coulomb gauge coincides with the 
real part of the effective potential governing the large-time evolution of the pair propagator. On the 
other hand, by looking only at the imaginary-time correlator at T = /3 , one looses any information 
about the collisional broadening suffered by the pair. 



4. Conclusions 

We addressed some very general aspects related to the propagation of a heavy QQ pair (used 
as an external probe to study medium properties) in a hot plasma. For the sake of simplicity we 
focused on the QED case. Medium corrections to the propagation of long wave-length modes are 
correctly described by the HTL effective action, which, for a QED plasma, turns out to be gaussian. 
This allows to get, for a static heavy-quark pair, an exact exponentiation of the resummed photon 
propagator leading to a closed Schrodinger equation governing the temporal evolution of the QQ 
propagator, with an effective potential containing both a real (related to screening) and an imaginary 
part (related to soft collisions, i.e. Landau damping). 
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At least two questions require further investigation. The first one is to study whether an 



exponentiation similar to the one we found holds also in the QCD case [15]. The second one is to 
consider the propagation of a finite-mass QQ pair, checking whether the ansatz of a closed 
Schrodinger equation for G > (t) is justified or a more refined strategy is required. This is presently 
the object of our efforts. 
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